Berry's Phases of Ground States of Interacting Spin-One Bosons: 
Chains of Monopoles and Monosegments 
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We study Berry's connection potentials of many-body ground states of spin-one bosons with 
antiferromagnetic interactions in adiabatically varying magnetic fields. We find that Berry's 
connection potentials are generally determined by, instead of usual singular monopoles, linearly 
positioned monosegments each of which carries one unit of topological charge; in the absence of 
a magnetic field gradient this distribution of monosegments becomes a linear chain of monopoles. 
Consequently, Berry's phases consist of a series of step functions of magnetic fields; a magnetic 
field gradient causes rounding of these step-functions. We also calculate Berry's connection fields, 
profiles of monosegments and show that the total topological charge is conserved in a parameter 
space. 
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I. INTRODUCTION 

In the presence of time-dependent periodical poten- 
tials V{t) = V{t -\- T) with a period T, an eigenstate 
of the instantaneous hamiltonian acquires a geometric 
phase known as the Berry's phase0. The Berry's phase 
is a global characterization of an eigen vector when it is 
transported in a Hilbert space. For a period of T dur- 
ing which an eigen state '^{y{t)) evolves adiabatically, 

the Berry's phase is simply <i>s — —Xm ( vjf \dt\ '$)dt. 
From the point of view of fiber bundles, the Berry's phase 
can be considered as the holonomy of a Hermitian fiber 
bundle; its base space corresponds to a parameter space 
and a fiber is an eigenstate which defines a mapping from 
the parameter space to the Hilbert space. This point of 
view was illustrated in 0- Berry's phases have been ob- 
served in a variety of experiments such as NMR and ro- 
tation of Hght polarisation in optical fibres etc0,0,|^|3. 

In condensed matter systems where there are large 
numbers of particles interacting with each other one way 
or the other, the subject of Berry's phases, or more gen- 
eral geometric phases becomes more fascinating and in- 
triguing. There are at least two interesting aspects of this 
subject: i) The effect of geometric phases on correlations; 
ii) The effect of correlations on Berry's phases. 

The first issue was addressed in quite a few different 
occasions. The best known result perhaps is the geomet- 
ric phases' effect on spin correlations. As pointed out 
a while ago for antiferromagnets, a Berry's phase dis- 
tinguishes integer and half-integer spin chains, or even- 
integer, odd-integer and half- integer spin square lattices 
and results in different ground states and excitations 
0, S I3- The other examples perhaps are the pos- 
sible geometric phases' effects on statistical transmu- 



tation and fractionalization in spin correlated systems 
[Tol [m ll^ IT3II . In mesoscopic magnetic systems ef- 
fects of Berry's phases on collective quantum tunneling 
were investigated [Illll El EHIl IH ; quantum in- 
terferences of Berry's phases were observed in molecular 
magnets po|. 

In this article we are going to address the second as- 
pect, correlations' effect on Berry's phases. We study 
Berry's connection potentials (one- form). Berry's connec- 
tion fields (two-form) and topological charge densities for 
ground states of spin-one bosons with antiferromagnetic 
interactions. In a parameter space 



{X4,a= 1,2,...M, 



(1) 



the connection potential A^, two- form connection fields 
Fab, and topological charge-current densities p, J are de- 
fined as 



Aa({^J;*) 

F,fc({Xj;v|/) 
Anp — — 



^abc 



IT 1 

< 



dXb ' dXa 



>}, 
(2) 



Consider Berry's connection potentials, two-form con- 
nection fields and topological charges when N spin-one 
bosons are in a magnetic field B = {Bx,By,Bz). For 
ground states of N-noninteracting spin-one bosons, one 
easily confirms that the two-form connection fields are 
monopole-like due to an N-fold magnetic monopole at the 
point of B = in the parameter space of {Ba}, a = 1, 2, 3 
or X, y, z. And the topological current is zero. That is 



'Permanent address 



ba = leabcFbc - Q(B)^,g(B) = N 
p = g(B)5(B),J = 0. 



(3) 
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We should emphasize here that the topological charge is 
localized at the origin B = in the parameter space. 
Note that the Berry's two-form field b is a function of 
external magnetic fields B. 

For N spin-one bosons with antiferromagnetic interac- 
tions, we are going to show that the profile of topolog- 
ical charge density is a linear chain of monopolcs. In 
the presence of a field gradient, we also find nonsingu- 
lar monosegments where topological charges distribute 
smoothly along certain direction and topological cur- 
rent flows. And the density profile is a linear chain of 
monosegments. Details of topological charge density pro- 
file depend on the even-odd parity of N, magnetic field 
gradient and spin relaxation. 

The plan of this article is as follows. In section II, we 
introduce the system which interests us and describe the 
algebras needed for this investigation. In section III, we 
calculate the Berry's connection fields of N interacting 
spin-one bosons in a uniform magnetic field. We show 
that antiferromagnetic interactions in general expel the 
topological charges outward from the origin and result in 
a linear chain distribution of monopolcs in the parameter 
space. In section IV, we study the effects of the spin non- 
conserving process; we calculate the monosegment profile 
in the presence of a magnetic field gradient. In addition, 
we address the Landau-Zener effect. 

In both section III and IV, we show that Berry's phases 
are suppressed because of antiferromagnetic interactions 
and as functions of magnetic field consist of a series of 
step-like functions. Furthermore, the shape of steps is 
determined by a magnetic field gradient. We note that in 
the absence of interactions, the many-body Berry's phase 
is simply the sum of each individual spin-one bosons and 
doesn't depend on magnetic fields. 

In section V we outline an alternative description of 
spin-one bosons using quantum rotor models. In the ef- 
fective representation, we show that the problem of in- 
teracting spin-one bosons in a magnetic field gradient is 
equivalent to a quantum rotor coupled to a quadrupole 
field. 



In a single mode approximation, the creation and an- 
nihilation operators are defined as 



(5) 



a = m = 0, ±1. Xo(r) is the lowest orbital mode and 
I 1,to) is a spin-one state with — m; S^fs are three 
matrices; 




= _ 1 1 sy 



5" 




(6) 



The dilute limit of spin-one bosons with antiferromag- 
netic interactions was studied in a few experiments in 
pH 12^ and also theoretically investigated in various 
works |23. 13 HE US' A geometric-description-based 
nonperturbative approach to strong-coupling limits was 
proposed in [23, Hg. In optical lattices, correlated Mott 
states of spin-one bosons were studied in a series of re- 
cent papers |2^ HD, |^ 113; ^P™ nematic, spin singlet 
and dimerized-valence-bond crystals were found for high 
dimensional and low dimensional optical lattices. Uncon- 
ventional spin disordered condensates in homogeneous 
limits were proposed in various papers [2^ H^; many 
properties of spin singlet condensates were further ex- 
plored in low dimensional optical latticesjs^. In fast 
rotating traps, correlated quantum liquids of spin-one 
bosons were also investigated jsE 113 . 

It is convenient to introduce the following creation- 
annihilation operators 



^1 



1 

i 
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{^-1 + 4) 



(7a) 

(7b) 
(7c) 



II. THE MICROSCOPIC HAMILTONIAN 



In this representation total spin operators are defined as 



We consider spin-one bosons in an optical trap in the 
dilute limit defined by na^ ^1 where a is the s-wave 
scattering length and n the density. The hamiltonian is 
given as 

H = I 1^ V^t (r) . (r) + C/(r) (r)^„ (r) (4) 
+ 7B .^t(r)s„^^^(r)+ |v3t(r)^t,(j.)^3^,(j.)^^(r) 
+ y Vii(r)^l,(r) S^^ • S„,^, #,(r)^0(r)} . 

Interactions between atoms are approximated by spin- 
dependent contact interactions. 



S"=Vits^^^^^S^^ = -ze„;3^; (8) 

"0^, V'a (from now on, a = x,y,z) are usual bosonic 
operators obeying the following commutation relations: 

Using these results, the hamiltonian can be written in 
the following form 

H = iea-9o-2g2)N + goN^ (10) 

+ 92 Stot + 7B • Stot 
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with 



90 = Y J'^^\^o{r)\^ 



(11) 



92 



C2 



dr|xo(r)|^. 



In terms of singlet pair creation operators A, 

it|0) = -^(-2^1iVl+V'5V'j)|0), 

the total spin operator can also be written as 
SL = ^(iV + l)-6iti, 



(12) 



(13) 



with N = 'tp\xpoi the number operator. Consequently, we 
rewrite the Hamiltonian as 

H = {eo-go-g2)N + {go + g2)N^ (14) 
- 632(i^i)+7B- 

A. Spectrum 

The Hamiltonian commutes with the number opera- 
tor, squared magnitude of the total spin operator S^^j 
and operator A^A. Eigenstates of the hamiltonia-n are 
simultaneous eigenstates of the number operator A'^ and 
operator A^A. If \N,Ns) is a state with a total number 
of A'' particles and Ng pairs of singlets, then we have 

SLi N,N,) = {N- 2N,){N - 2N, + 1)| N, ) (15) 

ov S = N -2Ns. 

For N particles in a magnetic field B along ^-axis, the 
spin Hamiltonian is 



H = g^Sl, + jB-St 



(16) 



Obviously, an external magnetic field along the z-axis 
splits the 25+1-fold degeneracy of states with given spin 
S and remains to be a good quantum number. The 
energy of an eigen state \S, Sz = m > is 



Es,m = g2S{S+l)+-fBm. 



(17) 



The lowest energy state for a given spin S is \S,—S). 
The low energy spectrum at different magnetic fields can 
be found in FIG.l. 
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FIG. 1: Energy spectrum as a function of magnetic field. £s,m 
is given in units of 32, B in units of 52 7~^. a) is for an odd 
A'^ and b) is for an even N. 



For the following values of magnetic fields 



Bk = ^ (4fc + 5) 

7 



(19) 



with fc = 0, 1, 2, . . . , Mjv where Mn = ^ for odd N. 
we have levelcrossings in groundstates. 



N even 



B. Levelcrossings 

1. N odd 

Levels 1 51,7711 > and \S2,m2 > cross at 



B = 



52 52(52 + l)-5i(5i + l) 



7 



mi — 1712 



(18) 



For an even N, the levelcrossings in the groundstate 
take place at 



Bfe = ^(4fc + 3), 
7 



(20) 



with fc = 0, 1, 2, . . . , Mn where Mn = for even N. 
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III. MONOPOLES AND BERRY'S 
CONNECTION FIELDS OF N INTERACTING 
SPIN-1 BOSONS IN HOMOGENEOUS 
MAGNETIC FIELDS 



In a spherical coordinate system, wc denote a magnetic 
field B, of magnitude B pointing in the direction of n = 
(cos (j) sin 6, sin ^ sin 6, cos 6) as 



A. Local connection fields and topological charge 
densities 



(23) 



The microscopic many-body state oi \S,Sz = —S) is 
given as 



-S) = C(Vili)^i^lO). 



Here A is the singlet creation operator defined before; C 
is a normalization factor 



with 



C = 6- 



1 



/(M, S) =S\M\2 



m{ 2M + 2S+1)\\ 
(25+1)!! ■ 



(21) 



For a magnetic field B pointing at n direction, ground 

states 



\9 >= *s(s),-s(B)(n) 



(24) 



are states \ S,Sz = —S;n) with the ^-axis defined along 
n-direction. The corresponding many-body wave func- 
(22) tions are 

I 



*s(s),-S(B)(n) = C 



n S 



N-S , 

A^ — \0). 



r 



(25) 



And finally S* is a function of _B, the magnitude of 
external magnetic fields. As shown in FIG.l, 



Ml-. 



S = So{N) + 2 J2 ©(|B|-Bfc) 



(26) 



fe=0.1,2.. 



Sq{N) is zero for an even N but is unity for an odd N. 

For these correlated ground states, Berry's connection 
potentials in spherical coordinates {B,9,(j)) can be de- 
fined as 

d 

A(B) = -Jm < *s,-s(n)|^|*s.-s(n) > es 

1 d 
-^1m < *s,_s(n)| — |4's,_s(n) > eg 

-sLe^'' < *5,-5(n)|^|*s,-5(n) > e,. (27) 
A direct calculation of A yields desired results 

cosgQ(|B|) 
sine |B| 

Mm 

Q{\B\) = qo{N) + 2 ^ Q{\B\-Bk) (28) 



qaiN) = 1 for an odd number of particles and q{](N) = 
for an even number of particles. The two-form connection 
fields are 



Q(|B|)- 



B|2 



1 9Q(|B|) cos6i 
'B dB sva.6 



(29) 



It is worth emphasizing that the radial component of 
b-fields can be attributed to multiple shells, each of which 
carries exactly two units of charges and is located at B = 
Bk- Indeed, if one defines Anps = V • (b • 63)63, one 
obtains 



PB = qo{N)SiB) 



1 

2^ 



(30) 



k=0,l,2.. 



/c=0,l,2... 



which indicates an isotropic topological charge distribu- 
tion of a series of shells ai B = B^. 

The topological charge due to the ^-component of b- 
fields on the other hand consists of two contributions: 
a) an isotropic charge distribution of spherical shells ex- 
actly identical to —ps given above; b) linearly distributed 
monopolcs located at B = zLBke^ each of which also car- 
ries one unit charge. As a result, the total amount of 
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topological charge due to 6'-coniponent vanishes identi- Finally, the total topological charge and current den- 
cally on each shell; however, the ^-component deforms sities are 
the isotropic distribution on each shell completely. 



qoSiB) 



1 



27rsin6l|B|2 



J2 S{\B\ ~ BkMe ~ n) + 5{e)] 



fc=0,l, 



1 d^Qi\B\) 
47r|B| dB^ 



cot{6)e^ 



(31) 



Compared with noninteracting systems, we find follow- 
ing three new features in Berry's connection fields and 
topological charges. 

i) Following Eas l29l3Tl we find that the total topo- 
logical charge within the kth shell bounded by an outer 
surface Bi — Bk + h and an inner surface B2 — Bk — h 
where Bk+i — Bk ^ h is a. conserved integer independent 
of B and index k. Indeed, 

— / B • egdS - — / B • ends = 2. (32) 

Here the integration is over two surfaces 5i.2 of _B = Bi 2- 
Therefore each shell defined above still carries over all 
two units of topological charges. If we introduce the total 
charge of a surface with radius B as 

<t{B) = 2vry p{B,9)B^ sinO dO, (33) 

we find 

a{B)^qo{N)SiB) + 2 ^ 5(|B|-Bfc). (34) 

fc=0,l,... 

At last, topological currents on each shell circulate 
around the z-axis. 

As shown in the introduction, for noninteracting spin- 
one bosons, the topological charge is located at the center 
of the parameter space. Following Eas l31l3ll it is ev- 
ident that antiferromagnetic interactions between spin- 
one bosons lead to expulsion of topological charges from 
the origin B = 0. As a result, topological charges dis- 
tribute on different shells with radius Bk, k — 0, 1, ...Mn\ 
each shell carries two units of charges. 

However, following Eg 1311 one also confirms that the 
total topological charge Qt is conserved; Qt is exactly 
the number of particles in the many-body state, indepen- 
dent of interaction strength. That is 

Qt = / p{B)dB = N. (35) 



One can easily show that the solutions in Eas l34l351 are 
independent of the specific class of ground state wave 
functions chosen for this investigation, though the results 
in Eg 13 II do depend on choices of wave functions. 

ii) Furthermore, on each shell the distribution of b- 
fields breaks the rotational symmetry, unlike in non- 
interacting cases; consequently the charge distribution on 
each shell is highly anisotropic. The expulsion of topo- 
logical charges for the particular set of states studied here 
is along the z-axis only. One can verify that in the case 
of homogeneous magnetic fields two units of charges on 
the fcth shell are located at 9 = 0,71, B = Bk points; that 
is, all charges on each shell are carried by two monopoles 
located at the northern and southern poles of shells. So 
the density profile is a chain of monopoles located at 
B = ±Bfce^, fc -I- 0,1,. .Mat. 

This structure however is not generic. As we will show 
in the next section in the presence of a field gradient, 
charges on each shell are carried by smooth structures of 
monosegments instead of monopoles. 

iii) On each shell oi B = Bk, k — 0, l,2...Afjv as 
shown in FIG. 2, the b-field has a new component along 
^-direction. This is due to the level crossing or more pre- 
cisely presence of different spin states in exact ground 
states. The new component represents a Berry's phase 
when a closed path in the parameter space crosses B = 
Bk surfaces. We will come back to this point in the next 
subsection. 

B. Berry's phases I: rotating fields 

When a magnetic field B with a given magnitude ro- 
tates around the z-axis, a correlated state of spin-one 
atoms evolves along a path C shown in FIG. 3a). Dur- 
ing each period the state acquires a many-body Berry's 
phase. Consider a rotating magnetic B(i) = Bn{t) with 
n{t) — {cos (f){t) sin 9, sin (t){t) sin 9, cos 6); 6* is a constant 
and (f> — 2'KD,t is time-dependent. 

Previous results on connection fields indicate that 
Berry's phases of many-body states 

$b(C) ^ j K-dB (36) 
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By 



i 1 i 



- • - 

k f ' 



' r , , ^ 



-1 



b) 



r 1 . 
1 



-1 



c) 



FIG. 2: Schematic of topological field and charge distribution 
in the parameter space {Bx, By, Bz)- In a) we show the topo- 
logical field and charge distribution for Ai'-noninteracting spin- 
one bosons. In b) and c) we show the field and charge distribu- 
tion for an odd and even number of spin-one bosons with an- 
tiferromagnetic interactions respectively; chains of monopoles 
are shown here explicitly. The b-fields only have radial com- 
ponent except on shells of B = Bk, k = 0, l,2...Mjv. Note 
that monopoles in b) and c) are also distributed on a series 
of shells with radius B = Bk\ and for the even case, there are 
no monopoles at the origin of the parameter space. 




By 



Bz 



By 



FIG. 3: Paths in a parameter space, a) A path corresponding 
to a rotating field; b) A path for a modulating field discussed 
in this subsection; the path is in 6 = 6o surface. 



depend on the magnitude of magnetic fields. Following 
calculations in Appendix C, 

$B -27r cos6iQ(B). (37) 

As shown in Eg 1371 at a given magnetic field which is 
much smaller than JBmjv> the Berry's phase is always 
strongly suppressed because of antiferromagnetic corre- 
lations. We summarize our results in FIG. 4. 

The non-analytical behavior of Berry's phases in rotat- 
ing fields is consistent with an anisotropic distribution of 
monopoles discussed in the last section. Consider two in- 
finitesimal paths Ci_2 centered at and also oriented along 
the z-axis. Ci is slightly above —BkGz and C2 slightly 
below. The Berry's phases evaluated in this way are sin- 
gular at points —Bke^ and experience a jump 



$b(C2) -$b(Ci) =4^. 



(38) 



Thus the step-function plotted in FIG. 4 indeed implies 
that the surface bounded by path Ci and C2 enclose a 
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FIG. 4: A Berry's phase for a rotating field witli — ^ and 



various B; B is given in 327 
number. 



N is taken to be an odd 



monopole at B = —BkGz- One can apply similar argu- 
ment at B = -BfeGx and arrive at the same conclusions. 



C. Berry's phases II: modulating fields 

For non-interacting cases, b-fields only have e^- 
components. If a path completely lies in a plane parallel 
to e^, a quantum state doesn't acquire a Berry's phase 
because the topological flux threading the path should 
be zero for an obvious reason. 

One of interesting aspects of Eq|251 is that there is a 
0-component in b-fields because of antiferromagnetic in- 
teractions. This yields a new possibility to study Berry's 
phases which is uniquely associated with interacting par- 
ticles. Consider a path that lies in a plane oi 9 = 9q, and 
is bounded by 



as shown in FIG. 3b). 

Following Eq[551 one obtains 



(39) 



' COSf^o [<p2 



i)[Q{B2) - Q{B,)]. 



(40) 



This Berry's phase is nonzero only when B = Bi,2 sur- 
faces are at two sides oi B — Bk surfaces. For B^+i > B2 
> Bk> Bi> Bk-i and 6*0 = 37r/4, 0b = (0i - 02)V2. 
The Berry's phase for modulating fields of this kind van- 
ishes identically for noninteracting particles. 

Discussions on geometrical phases are valid when slow 
spin relaxation is allowed so that the system can al- 
ways reach true ground states in different magnetic fields 
within practically relevant time intervals. Furthermore 
we also assume that the quantum symmetry restoring 
time is much shorter than measurement time; the issue 
of symmetry restoring of condensates of spin-one bosons 
was addressed in some detail in Ref.jSJ]. 



IV. MONOSEGMENTS AND BERRY'S PHASES 
OF N INTERACTING SPIN-1 BOSONS IN 
INHOMOGENEOUS MAGNETIC FIELDS 

A. Spin conservation and magnetic field gradient 

In a homogeneous magnetic field, because of total spin 
conservation there is no mixing between states with dif- 
ferent spins. However, a gradient in magnetic fields, as 
we will see doesn't conserve total spin and does mix states 
with different spins. This results in "level repulsion" 
when two spin states approach each other. 

To show that an inhomogeneous magnetic field violates 
the conservation of total spin, we consider the commuta- 
tor [iJ, Sfgf\ with H given by 

H .92 ^tot +7 JdiiB{ii)- 



Sa _ -a 

One easily verifies that 



■'tot 



dxB(x)-V'lS„^V'/3,S2 

V'a(y)Saa''0a'(y)) 

2 /dx7/.t(x)B(x) •S„0^^(x). 



(41) 



(42) 



This commutator only vanishes in a homogeneous mag- 
netic field but is nonzero when a field gradient is 
present(see Appendix B). 



B. Mixing of many-body states with different total 

spins 

For spatially varying magnetic fields B(x), we define 
new creation and annihilation operators ip'^, ip'^^ in a 
local triad where the local magnetic field points at the 
z direction; the local field operators can be obtained by 
the following spin rotation 



^aW = C/„^(x)?^^(x), 

with C/jj^(x) given by 

Ua/si-x.) = exp(i0(x) • Sq/j). 



(43) 



(44) 



The Zeeman splitting in terms of spin-rotated opera- 
tors is: 



Hze 



7 /dxS(x)^^t5^^^^^^.^^„,,.^. (45) 



with -B(x) — |B(x)|. On rotated spin-one fields, external 
magnetic fields always act along Z-axis. 
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One also finds that the spin-dependent and spin- 
independent interaction terms are invariant under the lo- 
cal spin rotation. However, the kinetic energy transforms 
nontrivially and acquires a new term, Hkin Hkin + Hi 
(see Appendix C). 

Consider a field distribution 



B(x) ^ Bo{l- G'z)z + G'Boxx. 



(46) 



We assume that G' is very small compared to the dimen- 
sion of the system, i.e. GTia <C 1. This ansatz was 
employed for the study of field gradient effects "2^. 0, 
which is determined by the relative orientation of mag- 
netic fields with respect to the z-axis, can be calculated 



states. We will study the resultant spectrum in the next 
subsection. 

C. An effective Hamiltonian close to crossing 
points Bk 

For values of |B| close to the level crossing points Bk, 
we have a small tunneling term calculated above. If two 
nearly degenerate states are far away from other lev- 
els, the Hilbert space can be truncated into a two-level 
Hilbert space and the effective Hamiltonian in this sub- 
space can be easily determined. 

When B is close to Bk, we therefore obtain an effective 
Hamiltonian in the truncated space 



6»(x) = z X B(x) = G'xy. 



(47) 



Following Appendix C, we find that the spin rotation 
considered above effectively orients external fields along 
the z-axis and results in a new term Hi in the single 
mode hamiltonian, i.e. 



and 



H[B{x)]^ H[B{x)e,]+Hi 



(48) 



(49) 



where e 



The tunneling matrix elements can be calculated as 



;i,-i|i7i|3,-3; 



N- 



N (50) 



For general level crossings in ground states, we obtain 
the following matrix elements. For an odd N, 

Afe = ( (2fc + 1), -{2k + 1) \Hi\ {2k + 3), -{2k + 3) ) (51) 

_ h^G'^ I {2k + 2){N -2k- \){2k + 3)(iV + 2k + 4.) 
Am y 

And for an even N, 



{Ak + 5){4:k + 7) 



Ak ^ {2k,-2k\Hi\2k + 2,-2k-2) (52) 
h^G'^ I {2k + 1){N - 2k){2k + 2){N + 2k + i) 



4m 



{Ak + 2,){Ak + b) 



Therefore a field gradient in magnetic fields has non- 
vanishing matrix elements between states with different 
spins and leads to mixing of corresponding many-body 



H = 



Es,^^s,{B) Ak 

Afc Esi^+2-Sk-2{B) 



(53) 



Sk — 2k + l for an odd N and Sk = 2k for an even A^.And 
once again, fc = 0, 1, 2. ..Mm- We have assumed that the 
gradient is small and ^ 52 so that the truncation 
discussed here is always applicable when B is close to 
Bk- 

For two levels |1, —1 > and |3, —3 > ai B ^ Bq 



092 

7 ' 



H 



Ei^^i{B) Ao 

^0 -E'3,-3(^) 



with 



Ao 



3 h^G'- 



70 2ra 



N- 



(54) 



(55) 



i^i.-i = £"3,-3 = Eo with Eo = —Sg2 when B = Bq- 
Introducc SB = B — Bq, &Ei,^i = i?o. One obtains 

SEi^-i = —jSB, 5i?3,_3 — —3jSB. The eigen values can 
be expressed as 



E±{6B) = -3.g2 - 2jSB±{-/'^SB^ + Aly- 



(56) 



D. connection fields of spin mixed many-body 
states 



Because of the field gradient, matrix elements of 
the Hamiltonian between states of different spins are 
nonzero. Eigenstates are generally superpositions of 
states with different spins. 

Close to level crossing points Bk, in the truncated 
Hilbert space the eigen states are 



*(n) - Ski I Sk, -Sk;n) + Sk2 \ Sk + 2, -Sk - 2; n). 



\Ski 



\Sk2\'^l; 



(57) 



the corresponding many-body microscopic wave func- 
tions are 
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^in)=SkiC{rl,f'' (it)^|0) + Sk2D{rlif''+^ 



where C and D are normalization constants. is a 

creation operator defined in a local frame discussed in 
Appendix CI. For a given k, 



Sk = 2k + 1 for an odd N 
Sk = "ik for an even N 



(59) 



Finally, as the magnetic field increases from below Bk 
to above Bk , Ski varies from one to zero and Sk2 from zero 
to one. Taking into account the effective Hamiltonian in 
the truncated Hilbert space derived in section IV C, we 
obtain the following field dependence of the coefficients 
<5/ci,fe2 when B is close to Bk 



Ski 



Sk2 = 



i^SBk - ^Ai+j^SBl) 



2{^HBl + A2 - ^SBk^^i+i^SBl) 



(60) 



2{^HBl + AI- ^SBky/Al+^^SBl) 



In Eq|^ we have introduced SBk = B — Bk- Ak {k — 
0, 1, 2...Mn) represent the matrix elements calculated in 
Eq.51; Ak = ((2fc+l),-(2fc+l) |iJi| (2fc + 3), -(2fc + 3) ) 
for an odd TV and = ((2fc),-(2fc) \Hi\ (2fc + 2),-(2fc+ 
2) ) for an even N. 

Discussions on connection fields and topological charge 
densities can be carried out, similar to those in the pre- 
vious section. After some straight forward calculations, 
we find: 



cos I 



'Qs(|B| 



IBI 



(61) 



a(B) 

10[ 



(it)- 



(58) 




10 



15 



20 



25 B 



b) 



FIG. 5: The surface charge cr{B) as a function of magnetic 
field in the presence of a field gradient G' = 1.32cm~^. B 
is measured in units of (727"^. a) is for an odd number of 
particles and b) for an even number of particles. 



Assuming the field gradient is small, we find that Qg(|B|) 
is identical to (5(|B|) defined in section III A when B is 
far away from degeneracy points Bk] close to Bk when 
j\B — Bk\ are comparable to however. 



Qg{\B\) = Sk + 2\S, 



k2\ 



(62) 



which varies smoothly from Sk to Sk + 2. We want to 
emphasize that both Ski and Sk2 are real functions of |B| 
and independent of 9 and (f>. 

Following discussions in the previous sections, one ob- 
tains the two-form connection fields Fat or b 



Qg(|B |) 
IBI 



1 dQg{B) COS0 
B dB sin 9 * 



(63) 



The topological charge density is 



p-go(iV)^(B) 



1 



27r|B|2sin6i dB 
Correspondingly, the surface charge becomes 



[S{9 - n) + S{9)] 
(64) 



a{B) = qo{N)S{B) 



dB 



(65) 



which is analytical at Bk, k = 0,1, ...Mm because of level 
repulsion. Following Eas l61l621 far away from Bk, cr{B) 
is vanishingiy small. cr(-B) as a function of B is numeri- 
cally plotted in FIG. 5. 

Eas l63l651 again clearly indicate shell structures. In- 
deed, when the width of each shell 2h defined before 
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Eq|2|is much larger than a characteristic width Wk 
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(66) 



the charge Qfc enclosed in the kth shell centered at B 
Bk is 



(67) 



Qk approaches two units when the ratio between 2h and 
Wk becomes infinity. And only on these shells, b fields 
have a ^-component. 

Therefore, Eg 1641 shows that on each shell the charge 
distribution is highly anisotropic. All charges on each 
shell in this case are carried by two monosegments. Each 
monosegment carries one unit of charge and is centered at 
B = zbSfce^; a monosegment is smooth along z-axis, the 
typical width of a monosegment located at B = iSfee^ is 
approximately Wk ■ Topological charges overall distribute 
in a chain of monosegments instead of monopoles. Fi- 
nally, in the presence of field gradients, results in EalSSI 
are still valid. 

Before we leave this subsection, we emphasize that b- 
fields are gauge invariant under a usual U{1) gauge trans- 
formation. Though this is hardly surprising by virtual of 
the two-form construction, in appendix E we neverthe- 
less present an explicit calculation to illustrate this point 
in terms of many-body wave functions. 




FIG. 6; Eigenvalues as functions of magnetic fields close to 
Bo; G" = 1.32 cm-\ 



with V as the rate. In FIG. 6, we show the energy of these 
two levels as a function of magnetic field in the vicinity 
of Bq. The steps in the Berry's phases as a function of 
magnetic field become rounded because of level repulsion. 



E. Landau-Zener effects 

B{t) varies adiabatically. When i — > — oo, the ground 
state is | 1, — 1 ) and when t +oo the ground state is 
I 3,— 3). Since B changes adiabatically, for the most of 
time the system remains in the ground state; however, 
as the change rate is finite, the system also makes tran- 
sitions to an excited state | 1, — 1 ) at i ^ oo with small 
probability. 

Denoting the excited state by + and the groundstate 
by — , the transition probability can be calculated with 
the following formulaji 



V. LARGE N LIMIT 
A. An effective hamiltonian in large N limit 

As proposed in a few previous works, the problem 
of interacting spin-one bosons can be mapped into a 
constraine d ci uantum rotor model in the large N limit 
I^^HEIMll^Ull' And any microsopic many-body state 
can be expressed in terms of a wave function i^in) of a 
quantum rotor characterized by its direction n. In fact, 
an arbitrary wave function V'(ii) represents the following 
microscopic wave function of a correlated state 



W+- 



exp 



~2Xm / dT{E+{T) ~E_{t)) 
Jc 



(68) 



Here tq is the complex value at which E+(t) = E-{t) 
and C a curve from t —oo to t ^ +oo passing above 
tq. Substituting the results found in Eq.59 into the above 
expression, we get 



W+.^ 



exp 



exp 



1 Tm 



dr WA5 -I- 72u2t-2 



(69) 



7f 



dn7/'(n) |n >; |n >= 



iV + 1, 
2N\ 



(n^i^triO). (70) 



Therefore, a state where n is localized on a two sphere 
corresponds to a polar condensate; an S-wave of n rep- 
resents a rotationally invariant spin singlet ground state. 
More discussions about connections between the two rep- 
resentations can be found in js^ ] . 

The hamiltonian for spin-one bosons in |n >- repre- 
sentation is given by a quantum rotor model 



= 32 + 7 B • L. 



(71) 
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10 15 20 25 

a) 




10 15 

b) 

FIG. 7: Berry's phase in the presence of a field gradient G' — 
1.32cm~'^. B is measured in units of g2'y~^ ■ a) is for an odd 
A''; b) is for an even N. 



The total spin of spin-one particles Stot{= L) is a differ- 
ential operator 



d_ 
9n' 



(72) 



i.e. the total spin operator is the angular momentum 
operator defined on the two-sphere where n lives; L is 
also a conjugate variable of n. Wave functions further 
observe the following Ising symmetry 



(73) 



for any N; this property of many-body wave functions 
was identified and emphasized in previous works on ho- 
mogeneous gases of spin-one bosons |23, ^8] and on spin- 
one bosons in lattices Isil Is^ l33 | . 

A direct calculation in Appendix F indicates the fol- 
lowing relation in the |n >-representation, 



V'i^a/3V'/3 NUa ha/3 n/3 . 



(74) 



Therefore, the field gradient induced quadratic Zeeman 
term discussed in Appendix D results in the following 
coupling term in the quantum rotor representation (up 
to a shift) 



Hi 



(75) 



Here 



/i„0(r) = V0,.V%S^^S« ; 



and Qap is the nematic order parameter 



?Q/3 = n/3 — -5af3- 



(76) 



(77) 



For a magnetic field distribution given in this article, 
the /i-matrix is 



G'2(S^S^)„^; 



(78) 



therefore, the effective hamiltonian for the gradient term 
in the large N limit is 



2m 



(79) 



which was also proposed in an early paper |28| . 



B. Effects due to anisotropy 

Taking into account anisotropy, the quantum rotor ef- 
fective hamiltonian for a magnetic field distribution given 
in Eq.Jinil is 



H ^ g2L^ +jB-L + ha(3QapN, 



2m 



(80) 



The tunneling term between different total spin states 
is given by 

{sm\hai3iplTp(3\s' m') = hapTapism, s'm') N 
Taj}{sm,s'm') ^ JdnY*^{n)nanpYs'm'{n). (81) 

For states | 1, — 1 ) and | 3, ~3 ), we have 



T„/3(3,^3;l,-l)-^/- 



1 i 
i -1 



(82) 





And the matrix element between I 3, — 3 ) and 1 1, — 1 ) is 



Nha0Tap{3,-3;l,-l) = ^^\[^N (83) 

2m V 70 

Close to Bq, states 1 1, — 1 ) and | 3, — 3) are nearly de- 
generate as mentioned before. In the two-level subspace, 
the Hamiltonian is identical to that in Eq.53; Ao(= hafs 
Taf)) is calculated above. One can then study the eigen 
states in this subspace; after redefining Aq, one obtains 
results identical to those in section IV C. 
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C. Berry's phases 

In the quantum rotor representation, the spectra of 
the Hamiltonian subject to Ising symmetries in Eg 1731 
are identical to those discussed in section II. The Berry's 
phases of many-body ground states in a rotating field 
can be evaluated in a straightforward way. All results 
in section III, IV can be easily rederived; moreover, this 
effective description gives a simple geometric interpreta- 
tion of results obtained in a microscopic calculation. We 
do not present detailed calculations here. 

VI. CONCLUSION 



contracts OOCCSPPIO, 02SIC25 and NWO-MK "projec- 
truimte" 00PR1929. 

APPENDIX A: BASIC ALGEBRAS 

One can verify the following algebra: 

[S",^^] = le'^P-'i,^ (Ala) 

[S",4>1] = ie^/^-^V^t (Alb) 

[S",S^] ^ ie^'^'^S'^ (Ale) 

[/5,?/'„] - -via (Aid) 

[P.i^i] = (Ale) 



In this article we discuss the Berry's connection fields 
of many-body states of spin-one bosons with antiferro- 
magnetic interactions. We show that unlike noninter- 
acting systems. Berry's connection fields are determined 
by a linear chain of monopoles; more over in the pres- 
ence of field gradients each monopole becomes a linearly 
extended object which we call a monosegment. Antifer- 
romagnetic interactions appear to result in anisotropic 
expulsion of topological charges from the origin of pa- 
rameter space and suppression of Berry's phases in small 
rotating magnetic fields. 

The peculiar Berry's phases of many-body states might 
be observed by studying resonance transitions between 
a ground state and a collective excitation, similar to 
NMR experiments carried out earlier In an adiabat- 
ically rotating magnetic field, the resonance frequency 
between these states should be shifted because of geomet- 
ric phases. For instance, at B less than Bq the shift in 
the resonance frequency of transitions between an excited 
state |2fc -I- 1, —2k — 1 > and the ground state |1, —1 > 
for an odd number of particles is 



5lj = —AkncosOoQ. 



(84) 



Here fl is the rotating frequency and is the angle be- 
tween the magnetic field and rotation axis z (see section 
IIIA for the geometry). 
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(Alf) 



Furthermore, we define the singlet creation operator: 

A'' = -^V^tv^t with a = a;,y,z (A2) 

= -L(v;jv;t_2^i^^t)_ 



This operator satisfies 



-(3 + 2/V). 



(A3) 



APPENDIX B: FIELD GRADIENT AND SPIN 
CONSERVATION 

In this appendix, we prove that a field gradient does 
not conserve the total spin of spin-one bosons. The sec- 
ond quantised form of the total spin operator Sj^j is given 
by 

SL = 2 ydx^t(x)^„(x) (Bl) 
We consider the commutator 



SL,7 y'dxB(x).V3t(x)S„/5V'/3(x) 



(B2) 



Calculations yield 
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'tot 



dxB(x) • -0t Sq/3-0/3,S: 

dx B(x) • 'ijj'lSaffi'p , 2 Jdxipl{x)il!a{x 

rfxB(x) ■ IplSa/ltpp , Jdx Jdy ■tpl{x)'tl!j^{y)Saa' ■ Spi3'tp/3'{y)lpa' i^^) 



(B3) 



r 



Using 



and 



[A, BC] ^[A,B]C + B [A, C] (B4) 



V'o(y)>^(x) =<5„^<5(x-y) (B5) 



we find that the first term vanishes. The second term can 
be simplified using the relations above and in the end we 
find results in Eq031 Now take B(x) homogeneous. We 
find 



'tot 



(B6) 



y"dxB(x) •i^ts„/3#,s: 
X (y"dyV^l(y)Sca'V'a'(y))-2B- Jdx^i{x)Sa,pi>f,{x 



= {-2i)B ■ (Stot X Stot) - 2B • Stot 

= {-i)B,e,,k[Sl^,Sl 
= 0. 



^tot\ 



So for a homogeneous magnetic field total spin is con- 
served and an inhomogeneous magnetic field breaks con- 
servation of total spin. 



APPENDIX C: CALCULATIONS OF BERRY'S 
CONNECTION FIELDS 



1. General discussions: without field gradient 

Define a new coordinate frame (x',y',z') with the 
2:'-axis in the direction of the magnetic field n = 
(sin0 cos (/), sin^ sin0, COS0), and the x'-axis in the (x,y)- 
plane. The creation operator in the local frame is 



(CI) 



Therefore, (^'5^-5(11) | d'^s-s{'^) ) is given by 



( vl/s,_5(n) I d*5.-s(n) )=C^S{Q e"*^ sin2(-)^t + ^i4> eos2(|) | ) , 



The results for connection fields follow this identity. 

Now the Berry's phase for N interacting spin-one 
atoms with total spin S is given by 

(j^Berry = -Tm <j> (*s,-s(n) |d*s,-s(n)) 
= -2713 co8{e). 



2. General discussions: with a field gradient 



Now for B - Bi. 
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; ^s,-sin) I dM>s,-sin) ) = -i cose[c^ Sk Sl,{B)d(t> { lA^iji" _,^'^''A^'^\ ) 

= -{5l^{B) Su + 5l^{B){Sk + 2)) cos Bdcp 

I 



3. Berry's local connection fields and topological 
charge densities 

As mentioned in the introduction, Berry's connection 
fields of tlie ground state are given by 



A(B)=-Xm(g|A|5) 



(C2) 



For a state 15, —S* >, following discussions in Appendix 
CI, 2, Berry's connection potentials in spherical coordi- 
nates (0, 0, B) are 



A(B)--^cot0e^ 



(C3) 



Now the Berry's connection fields are given by: 

(C4) 

cot(6i)ee. (C5) 



b = Vb X A(B) 

Q{B) , 1 dQ{B) 



"^^ ' B dB 
For B 7^ 0, 7^ vr and 6* 7^ 0, we find 

Anp = Vb • b 



(C6) 



J_d_,j.2QiB)_ 1 d IdQjB) 

B^dB^ 52 ''^ BsmOde^B dB 
1 dQ{B) 1 dQ{B) 



dB B2 dB 



= 0. 



At B = or 6* = TT, 6* = and B ^ Bk points, the 
b-fields are singular. Calculations around those points 
lead to results in Eal^ 



APPENDIX D: QUADRATIC ZEEMAN EFFECTS 
DUE TO FIELD GRADIENTS 

We show explicitly for the invariance of a spin- 
dependent term in the interaction. 
Let Uaf3 — exp[i0 • Sq,^]. 

Hint = y j'dxipl'iljjjSaa' ■ Sf3f3' ^a'i'l^' (Dl) 

2 



f/I,^,(x)^;,[/t,,,(x)^^, 



Now using 

Uaoi' Sq/^' Uj^r — 

{Uaa' X ■ Sa'j3' Uj^i Uaa' Y ' ^a'l3' ^^'/Ji 

= (x' • Sq;3 , y' • Sq/3 , z' • Sq/j) (D2) 

where the primes correspond to the rotated coordinate 
system and the fact that the inner product is rotation- 
ally invariant we find that the spin-dependent interaction 
term is invariant under the local spin rotation. 
Now the kinetic term 

(;2 



H = l^x -|-7^t(x)v2^„(x) 
becomes in terms of V^'J^, and -0^ 

y'dx^'t{/,,(x)V2t/i(x)^^(x) 

- /rfx^^V'l(x)V^V'l(x) 

+ l^x ~|^V''l(x)(7«/3(VC/t^).V^;(x) 



(D3) 



(D4) 



dx 



2m 

Using the relation 
d 



V''l(x)(C/„^(x)V2[/^^(x))0;(x) 



dxi 
and 



exp 



i6/(x) -S =i • S) exp 



de 

dxi 



dxi 



i9S 



(D5) 



(D6) 



for the ansatz in Eg 1461 we have H = Hq + Hi , and 

Hi ^ ^ jdxi:'iU^p{iG'Sy)p^U\^^^'^ (D7) 

- ^ jdxi/iu^fA-G'\sy)\^ul^^'^. 



By keeping only the zero mode contributions, we find 
■^i{Sy)lpi^f,. (D8) 



2m 
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This is the effective quadratic coupUng in the presence of 
a field gradient which was previously obtained in Ref.|26j|. 



APPENDIX E: GAUGE INVARIANCE OF 
TOPOLOGICAL FIELDS 

1. General 



We define the following basis 



\S,-S;n) ^C{S) e*^^^'* 



it^io). 



(El) 



r 



Close to Bk points, ground states are in general super- section IV D, wc find 
positions of these states; using the ansatz introduced in 

I 



n) = Ski{B)\ Sk, -Sk;n) + Sk2{B) e^'^''(«^*^^)| Sk + 2, -Sk -2;n] 

I 



(E2) 



with Ski.k2, dk functions of {B,(f>,l 
'5fci.fc2 as given in section IV D; dk 
Therefore we find 



Calculations yield 



0. 



dP 1 dP 

1 dP 

+ Qg{B) cos 9)e4 



BsinO 



(E3) 



Note that Berry's connection potentials A(B) transform 
in an expected form 



A(B) ^ A(B) + V SP{B, e, (j)] 
under a gauge transformtion P ^ P + 6P. 



h(B) 



1 dQg{B) cos 9 
B dB sin 61 ' 



(E4) 



(E5) 



We find on the k-th shell two monopoles, one at the 
southern and one at the northern pole. The connection 
field b(B) does not depend on the choice of P{B, 9, (p) 
i.e. the gauge. 



2. Case P{B,e,4>) =0 
We find in this limit, 



(E6) 



b(B) 



QgjB) 
52 



1 dQg{B) cos 9 
B dB sin 9 



ee (E7) 



APPENDIX F: ESTIMATE FOR COLD SODIUM 
ATOMS 

Consider a trap with N = 10® atoms, a density n = 
lO^'^cm^^ and a2 ~ ao = 0.32 nm. We have gp = — ^ for 
^■^Na. The Zeeman effect is characterized by 7 



1 = gpfJ-B 



-3.31 



(Fl) 



Note that 7 is negative. Previous results can be applied 
provided we take for the magnetic field B = —Bn{9,(p). 
One also finds 



Ao = 2.18 • IQ-^"^ K G'^ , 



(F2) 



where we have the following units: G" in m ^ and A in 
K. Furthermore 

32 = 10-15 K (F3) 
The eigenvalues of H in the vicinity of Bq are 

(SB) = ~3g2~ 2jSB (F4) 



±7\/4.23 • 10-21G2 to4 (^'4 +{SB)^. 
The relative level repulsion is 4 G'^ with G' in 
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APPENDIX G: CALCULATION OF MATRIX 
ELEMENTS IN A LARGE-N LIMIT 

Now the general operator tpl^hajs^/j with hafi some ma- 
trix acts on the wave function as foUows: 



{ipi{n)\ip'lhai3ip/3\^2in)) 



dm /dn2^^(0|(niv37)^V'l/ia/3#(riIV'l)'^|0)V'i*(ni)V'2(n2) 



dni / dn2 



2NI 

N + 1 



(Gl) 



2NI 

dn N nahapnp ^pl (n)V'2 (n 



N^nih^pnlri{ni)M^2){0\{n\i^^f-'{nji;X-'\0) 
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